Calhoun 

Imti'uliOMl  Archive  ■of  dw  NavJl  PostguduiHe  School 


Calhoun:  The  NPS  Institutional  Archive 
DSpace  Repository 


Theses  and  Dissertations 


Thesis  and  Dissertation  Collection 


1986-12 

Development,  validation  and  use  of  a 
computer-controlled  system  for  the 
investigation  of  phase  and  amplitude  shaded 
acoustic  arrays 

Butler,  John  D. 


http://hdl.handle.net/10945/21791 


Downloaded  from  NPS  Archive:  Calhoun 


DUDLEY 

KNOX 

LIBRARY 


Calhoun  is  a  project  of  the  Dudley  Knox  Library  at  NPS,  furthering  the  precepts  and 
goals  of  open  government  and  government  transparency.  All  information  contained 
herein  has  been  approved  for  release  by  the  NPS  Public  Affairs  Officer. 

Dudley  Knox  Library  /  Naval  Postgraduate  School 
411  Dyer  Road  /  1  University  Circle 
Monterey,  California  USA  93943 


http  ://w  w  w.  nps.edu/lEbrary 


BUI5LT  - 

X\Fi'iVA.L*  POS  ’Atux..  '  .lOOL 
MONTEREY,  CALIi  ORE  lA  9?Ci3  "002 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 


DEVELOPMENT,  VALIDATION  AND  USE  OF  A 
COMPUTER-CONTROLLED  SYSTEM  FOR  THE  INVESTIGATION 

OF  PHASE  AND  AMPLITUDE  SHADED  ACOUSTIC  ARRAYS 

by 

John  Dale  Butler 

December  1986 

Co-Advisors : 

A.B.  Coppens 
G.B.  Netzorg 

Approved  for  public  release;  distribution  is  unlimited. 

T230157 


SKUAifv  CiASWi?AT;ON  (*>F  Yh  S  PAGE 


REPORT  DOCUMENTATION  PAGE 


la  REPORT  SECURITY  CLASSIFICATION 

Unclassified 


lb  RESTRICTIVE  MARKINGS 


2  a  SECUR'Tr  CLASSIFICATION  AUTHORITY 


2b  DECLASSiFiCAT  ON  /  DOWNGRADING  SCHEOULE 


3  DISTRIBUTION  /  AVAILABILITY  OF  REPORT 

Approved  for  public  release; 


distribution  is  unlimited 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


S  MONITORING  ORGANIZATION  REPORT  NuM8ER(S) 


6d  NAME  OF  PERFORMING  ORGANIZATION 

Naval  Postgraduate  School 


6b  OFFICE  SYMBOL 
(if  ipphabie) 

61 


)»  NAME  Of  MONITORING  ORGANIZATION 

Naval  Postgraduate  School 


Sc  ADDRESS  {Ci fy  State  *nd  /iPCode) 


7b  ADDRESS  {City  Stitt  a nd/IPCode) 


Monterey,  CA  93943-5000 


Monterey,  CA  93943-5000 


8a  NAME  OF  FUNDING  /  SPONSORING 
ORGANIZATION 


8b  OFFICE  SYMBOL 
(If  ipphciblt) 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


Be  ADDRESS  (Cry  S  fare.  and  /IP  Code) 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO 

NO 

NO 

ACCE  SS  GN  NO 

M  TjTlE  f include  Security  ClitSific ition) 

DEVELOPMENT ,  VALIDATION  AND  USE  OF  A  COMPUTEP-CONTFOLLFD  SYSTEM 
FOR  ^IIE  INVESTIGATION  OF  PHASE  AND  AMPLITUDE  SHADED  ACOUSTIC  APpa.^3 


:  personal  authoris) 


Butler.  John  D. 


•  )a  Type  Of  REPORT 

1  3b  T'ME  COVERED 

14  DATE  Of  REPORT  (Y**r  Month  Day) 

IS  PAGE  COuNT 

Master's  Thesis 

FROM  TO  

1986  December 

94 

6  Supplementary  notation 


’  T  COSATi  CODES 

f  ElD 

GROUP 

Su8  GROUP 

18  SUBJECT  terms  (Conf/nue  on  reverie  if  necessary  and  identify  by  block  number ) 

Computer-controlled  testing.  Phase  Shaded  Array,  Amplitudi 

Shaded  Array,  Beam  Pattern  Measuring,  Dipole  Radiation, 
Acoustic  Arrays 


9  ABSTRACT  (Conf/nue  on  reverie  if  neceiury  and  identify  by  block  number) 

A  computer-controlled  system  for  automated,  remote  measurement  of  phase  and 
amplitude  shaded  doublet,  triplet,  and  quadruplet  acoustic  arrays  beam  patters  was 
developed  and  tested.  Experimental  data  were  collected  with  an  anechoic  chamber  to 
minimize  backscatter  and  other  interference  effects. 


Preliminary  impedance,  reactance  and  resistance  plots  were  made  on  several 
speakers  to  ensure  that  a  matched  set  with  similar  characteristics  was  utilized  for 
the  arrays.  Once  the  arrays  were  established,  a  computer  program  was  developed  for 
an  HP-86  computer  which  allowed  it  to  rotate  the  arrays,  sample  and  collect  data,  and 
plot  the  resul ts . 


Various  combinations  of  amplitude  and  phase  shading  were  applied  to  the  dipole, 
tripole,  linear  quadrupole  and  quadratic  quadrupole  array.  The  results  were  then 


.0  O  S'R’SurON  ■  AVAILABILITY  OE  ABSTRACT 
KX.NCLASSifiEOlJNL  MlTED  □  SAME  AS  RPT 


□  otic  USERS 


21  ABSTRACT  security  ClASSif 'CATION 

Unclassified 


22a  NAME  Of  RESPONSIBLE  NOiViOual 

Alan  B.  Coppens 


22b  TELEPHONE  (include  Ares  Code) 

408-646-2941 


22c  OFF K  E  SYMBOL 

61Cz 


DD  FORM  1  473,84  mar 


83  APR  edition  may  be  used  until  exhausted 
All  other  ed't.ons  are  obsolete 


SECURITY  CLA SS.FiCATIQN  of  Th-S  PAG E_ 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PACE  (When  Data  Entered) 


19.  ABSTRACT  (continued) 

compared  to  computer  generated  beam  patterns  which  agreed  to 
uncertainty  was  mainly  generated  by  the  slight  deviations  in 
of  elements. 


±  2dB .  This 
the  positioning 


s  N  0102-  LF-  014-  6601 


Unclassified 


2 


SECURITY  CLASSIFICATION  OF  THIS  PAGEfWh#n  Dmf  Enfrmd) 


Approved  for  public  release;  distribution  is  unlimited. 


Development,  Validation  and  Use  of  a  Computer-Controlled 
System  for  the  Investigation  of  Phase  and  Amplitude 
Shaded  Acoustic  Arrays 


by 


John  D.  Butler 

Lieutenant  Commander,  United  States  Navy 
B.S.,  University  of  Texas  at  Austin,  1974 

Submitted  in  partial  fulfillment  of  the  requirements 

for  the  degree  of 

MASTER  OF  SCIENCE  IN  ENGINEERING  ACOUSTICS 

from  the 


NAVAL  POSTGRADUATE  SCHOOL 
December  1986 


ABSTRACT 


A  computer-controlled  system  for  automated,  remote 
measurement  of  phase  and  amplitude  shaded  doublet,  triplet, 
and  quadruplet  acoustic  arrays  beam  patterns  was  developed 
and  tested.  Experimental  data  were  collected  with  an 
anechoic  chamber  to  minimize  backscatter  and  other 
interference  effects. 

Preliminary  impedance,  reactance  and  resistance  plots 
were  made  on  several  speakers  to  ensure  that  a  matched  set 
with  similar  characteristics  was  utilized  for  the  arrays. 
Once  the  arrays  were  established,  a  computer  program  was 
developed  for  an  HP-86  computer  which  allowed  it  to  rotate 
the  arrays,  sample  and  collect  data,  and  plot  the  results. 

Various  combinations  of  amplitude  and  phase  shading  were 
applied  to  the  dipole,  tripole,  linear  quadrupole  and 
quadratic  quadrupole  array.  The  results  were  then  compared 
to  computer  generated  beam  patterns  which  agreed  to  ±  2dB. 
This  uncertainty  was  mainly  generated  by  the  slight 
deviations  in  the  positioning  of  elements. 
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I .  INTRODUCTION 


A.  PURPOSE  OF  STUDY 

The  purpose  of  this  investigation  was  to  develop  and 
verify  a  computer-controlled  experimental  system  to  measure 
and  plot  the  pressure  fields  of  various  acoustic  arrays  in 
an  anechoic  environment.  This  was  to  be  done  using  various 
combinations  of  phase  shading  and  amplitude  shading  of 
doublets,  triplets,  and  quadruplets.  Experimentation  was 
conducted  under  computer  control  to  maximize  the  amount  of 
data  available  and  to  ensure  the  highest  degree  of  accuracy. 
Additionally,  a  comparison  was  conducted  with  the 
theoretical  predictions  obtained  by  Park  [Ref.  1], 

B.  BACKGROUND 

The  use  of  dipoles  and  quadrupoles  for  beam  forming  and 
direction  finding,  both  in  air  and  in  water,  has  been  a 
subject  of  intense  investigation  since  the  1940's.  The 
basic  radio  direction  finding  apparatus  of  the  1940's  was 
nothing  more  than  a  set  of  orthogonal  dipoles  which 
consisted  of  two  crossed  loops.  The  relative  signal 
amplitude  from  the  two  dipoles  was  used  to  determine  the 
direction  of  arrival  of  a  signal.  The  angle  is  equal  to 

Arctan  (y/x)  where  x  was  the  amplitude  of  one  dipole  and  y, 

c 

the  amplitude  of  the  other  dipole.  This  presented  a  180 
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bearing  ambiguity  which  was  resolved  by  comparing  the  phase 
of  the  x  and  y  dipole  signals  with  the  phase  of  a  co-located 
omni-directional  antenna.  [Ref.  2] 

In  the  1950's,  extensive  studies  were  conducted  using 
planar  arrays  with  multiple  acoustic  and  electromagnetic 
dipole  elements.  These  studies  led  to  the  formulation  of 
planar  arrays  employed  in  radar  systems  of  the  1960's. 
Further  research  resulted  in  the  ability  to  steer  beam 
patterns  by  use  of  multiple  elements,  phase  shading,  and 
amplitude  shading. 

In  1976,  work  was  done  by  Hughes  and  Thompson  to 
investigate  acoustic  communications,  particularly  in 
underwater  acoustics  for  cases  where  the  position  of  the 
source  or  the  receiver  is  fixed.  They  investigated  the 
tilted  directional  response  patterns  below  the  horizontal 
plane  due  to  amplitude  weighting  and  a  single  90°  phase 
shift.  Here,  a  method  of  achieving  a  tilted  pattern  from  a 
symmetrically  shaped  array  without  using  numerous  phase 
shifts  or  time  delays  was  discussed.  [Ref.  3]  Moses  and 
Smith  also  examined  beamforming  this  time  with  a  cylindrical 
array  using  real  (rather  than  complex)  shading  coefficients 
[Ref.  4], 

In  1985,  Flanagan  examined  the  use  of  computer-steered 

microphone  arrays  which  were  two-dimensional  planar  arrays 

of  9  x  7  elements.  Here  an  array  was  designed  and  computer 

steered  to  discriminate  against  sound  arrivals  from  all  but 

the  desired  directions.  [Ref.  5] 
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One  example  of  the  Navy's  current  use  of  dipole 
technology  is  in  sonobuoys.  The  SSQ-53  DIFAR  Sonobuoy 
possesses  both  an  omnidirectional  and  bi-directional 
hydrophone  capability.  The  SSQ-77A  VLAD  Sonobuoy  possesses 
an  upper  and  lower  array  strung  in  a  single  line  array. 
Beamforming  is  accomplished  by  the  use  of  both  phase  shading 
and  amplitude  shading. 

These  advancements  in  the  study  of  arrays  ranging  from  a 
single  dipole  to  large  scale  planar  arrays,  coupled  with 
computer  control,  form  the  basis  for  this  further  study  of 
shaded  arrays.  There  is  currently  a  related  research  effort 
ongoing  in  the  study  of  computer  controlled  arrays  and  beam 
patterns  being  conducted  by  the  Environmental  Physics  group 
at  the  Naval  Postgraduate  School. 

C.  REPORT  OVERVIEW 

In  Chapter  II,  a  basic  introduction  to  theoretical 
pressure  and  directivity  calculations  of  a  dipole  is 
developed.  This  is  then  used  as  a  base  to  expand  to  the 
physical  case  and  the  electrical  linear  systems  case. 
Finally,  these  results  are  applied  to  different  forms  of 
quadruplets  and  triplets  to  show  theoretically,  the 
expected  results. 

Chapter  III  contains  a  summary  of  the  beam  pattern 
measuring  system,  the  computer  program  development, 
experimental  procedure,  and  representative  data  and  plots. 
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Additionally,  it  provides  a  description  of  the  system  which 


was  assembled 

for  measuring 

the  pressure 

fields  of 

the 

various  arrays. 

Chapter 

IV 

contains  a  summary  of  data 

analysis 

and 

ver if ication 

• 

Additionally , 

representative  plots 

are 

provided  for 

amplification  of 

arguments . 

Chapter 

V 

contains 

concluding 

remarks 

and 

recommendations  for  further  work. 

Appendix  A  provides  a  brief  description  of  each  piece  of 
equipment  and  summary  of  its  capabilities  as  applied  to  this 
experimental  investigation. 

Appendix  B  lists  the  BASIC  computer  language  program 
used  for  data  acquisition  and  formulation  in  this  study. 

Appendix  C  contains  specific  plots,  both  experimental 
and  theoretical,  obtained  as  a  result  of  this  investigation 
and  that  of  Park  [Ref.  1],  These  plots  were  obtained  for 
beam  patterns  of  various  array  configurations. 
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II.  THEORY 


The  dipole  has  been  extensively  explored  and  is  now 
utilized  as  a  basic  unit  in  the  construction  of  a  great  many 
arrays  and  antennae.  In  this  section  we  will  elucidate  the 
theory  using  two  approaches  1)  an  acoustic  approach  in  2 
dimensions  and  2)  a  linear  systems  approach  in  three 
dimensions.  These  two  approaches  are  related  to  each  other 
and  yield  the  same  result  after  minor  changes  are  made  in 
the  geometry  of  the  problem.  Both  are  presented  because  of 
direct  Navy  applications  dealing  deal  with  SONAR  arrays. 

Following  the  presentation  of  each  of  these,  the 
equations  for  the  triplet  and  quadruplet  arrays  utilized  in 
this  research  effort  will  be  shown. 

A.  DIPOLE  RADIATION-ACOUSTIC  APPROACH 

Because  it  is  usually  desirable  and  useful  to  be  able  to 
determine  the  direction  from  which  an  acoustic  signal 
originated,  arrays  were  developed  to  provide  this 
information;  a  simple  example  of  this  is  the  acoustic 
dipole.  A  simple  dipole  consists  of  two  simple  sources  of 
equal  strength  but  opposite  phase,  separated  by  a  distance  d 
Fig.  (2.1). 
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Figure  2.1  Geometry  used  in  Deriving  the 

Radiation  Characteristics  of  an 
Acoustic  Dipole. 


The  pressure,  p(r,0,t),  measured  at  the  field  point  (r,0) 

due  to  source  1  and  2  respectively  are 

Pl  =  |A/(r  +  Arj)]  exp[j(cot  -  k(r  +  Afj)]  ( i ) 

P2  =  [-A'(r  -  Ar2)]  exp[j(on  -  k(r  +  Ar2))  (2) 

where  A  is  the  pressure  amplitude  of  an  individual  source 
measured  at  one  meter  from  the  source,  r  is  the  distance 
from  the  field  point  to  the  midpoint  between  the  sources, 
Cl)=  2  77  f  is  the  angular  frequency  and  k  is  the  wave  number. 

The  minus  sign  associated  with  P2  accounts  for  the  180 

phase  difference.  The  total  acoustic  pressure  is  the  sum  of 
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Pi  and  P£ , 


P(r.0,t)r(A'r){([exp(-jkAr,)]  [ l+( Atj  r)|) -  {[exp(-jkAr2>]  [l»(Ar2/r)]]]exp(j(cot-  kr))  (3) 
Considering  only  field  measurements  where  the  separation  d  is 

small  compared  to  the  distance  r,  we  can  use  the  far  field 
approximations,  ^ri,2  =  (d/2)  sin  0  and  Arj  v'r  <<  i 

This  yields  the  acoustic  dipole  radiation  pressure  field, 

I’( r,G.t)  =  -j  (2A/r)  sin  {(kd  2)sin  0)}  exp  j(o)t  -  kr)  (4) 


(j) 


Figure  2.2  Cartesian  and  Polar  Representation 
of  Pressure  Amplitude  Distribution 
of  Dipole  Radiation.  [Ref.  6] 
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As  noted  in  Eq.  (4)  and  seen  in  Fig.  (2.2),  plotting  the  far 
field  for  d<<A,  pressure  nulls  exist  when  (kd/2 ) sin©  =77n 
where  n  is  an  integer.  [Ref.  6] 

The  twin  lobes  and  nulls  of  a  dipole  result  in  an  18fP 
ambiguity  when  trying  to  localize  a  source.  If  the  dipole 
is  combined  with  an  omnidirectional  receiver,  the  ambiguity 
can  be  resolved.  In  Fig.  (2.3),  various  combinations  of 
dipole  and  omnidirectional  receivers  are  combined  to  provide 
a  cardioid  response. 


Figure  2.3  Directional  Diagrams  of  Various  Combinations 
of  Bidirectional  and  Nondirectional 
Microphones.  [Ref.  7] 

Here  R,  is  the  voltage  output  of  the  omnidirectional 
receiver  and  R2  the  voltage  output  of  the  directional  dipole 
receiver.  It  is  not  important  for  the  two  sources  to  be  of 
the  same  sensitivity,  but  the  ratio  R^/Rx  must  be  unity  for 
true  cardioid  response.  The  phase  shift  between  these  two 
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receivers,  which  yields  directionality,  is  achieved  through 
a  difference  in  the  phase  characteristics  of  the  omni  which 
is  differentiated  (high  pass  filtered)  to  compensate  for  the 
linear  frequency  dependence  of  the  dipole  sensitivity.  [Ref. 
7]  . 

This  design  is  presently  used  in  SSQ-77A  VLAD  Sonobuoys 
where  a  cardioid  is  formed  pointing  in  the  downward 
direction.  This  maintains  the  null  in  the  upward  direction 
in  order  to  alleviate  high  frequency  interference  from  wind 
and  sea  surface  noise. 

B.  LINEAR  SYSTEMS 

Consider  the  geometry  for  a  dipole  as  shown  in  Fig. 
(2.4)  : 


Figure  2.4  Dipole  Geometry  [Ref.  8] 
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In  this  representation  as  developed  in  Ref.  (7) ,  Cj  and  C_x 
are  complex  weighted  point  sources  which  are  free  to  rotate 
in  the  x-y  plane  and  d  is  the  distance  between  elements. 

The  complex  weighted  response  factor  ci  t^ie  Point 

source  is 

ci  =  ^  exp(j<p.)  (5) 

where  ai  is  the  amplitude  and  ^i  P^ase  angle.  For  a 

dipole,  the  complex  weighted  response  factors  are  Cj  =  aj 

and  c-i  =  -3j  * 

For  a  wave  originating  from  a  distant  source  and 
propagating  to  the  receiver,  it  is  observed  that  there  is  an 
essentially  constant  pressure  magnitude  over  the  dipole 
aperture.  Because  the  phase  varies  with  position  along  the 
wave  front  in  a  traveling  wave,  the  response  of  the  dipole 
is  given  as  a  function  of  frequency  and  direction  cosine  u 
where  u  =  sinGcosNk.  The  response  is 

D(f,u)  =  Cj  exp(jTtdu'X)  +  c2  exp(-j7tdu,  X)  (6) 

This  equation  is  simplified  for  the  dipole  to: 

D(f,u)  =  2ja,  sin  (rcdu/X)  (7) 

Since  u  is  a  direction  cosine,  it  is  only  valid  over  the 
interval  -1<u<1.  As  a  result,  0<0<rc  and  0S\}/^2Jt  .  u  =  o 
implies  a  broadside  or  on-axis  geometry  and  u  =  1  implies 
an  end-fire  geometry.  A  plot  of  D(f,u)  as  a  function  of  u 
can  be  used  to  approximate  the  response  of  a  dipole.  [Ref. 
8] 
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An  example  of  this  two  element  response  of  two  point 
sources  with  a  separation  d  equal  to  one  half  the  wavelength 
of  the  traveling  wave  results  in  a  real  response  for  values 
of  u  from  (-X  2d)  to  ( +  X/2d)  •  If  the  sources  are  assigned 

dipole  weightings,  substitution  of  the  functions  ©  and  for 
direction  cosine  u  result  in  the  unnormalized,  far  field 
beam  pattern 

D(f,0.y)  =  j2a,  sin[(nd/X)  sin  0  cos  \\)}  (8) 

This  linear  systems  theory  can  be  easily  applied  to  the 
axial  quadrupole  and  Tesseral  quadrupole.  The  axial 
quadrupole  is  an  array  of  three  complex  weighted  point 

sources  where  c,  =  a,  =  c.t  =  a.,  =  1  and  cQ  =  aQ  =  -2  as  in 

Fig.  (2.5). 

Here  CQ  can  also  be  thought  of  as  2  sources  with  each  having 
an  amplitude  weight  equal  to  -1. 


Figure  2.5  Axial  Quadrupole  Geometry  [Ref.  8] 
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With  the  1,-2, 1  amplitude  weighting,  the  far-field, 

unnormalized  directivity  function  is  defined  as: 

D(f,fx)  =  [exp(-j27tfxd  2)]  -  2  +  [exp(j2rcfxd  2)]  (9) 

where  f'x  =  u,  X  =  sinGcosy  X  .  This  reduces  to: 

D(f.G.V(/)  =  2  [cos((md  X)  sinGcosy)  -  1]  (10) 

In  the  case  of  a  Tesseral  quadrupole  Fig.  (2.6),  a 
planar  array  of  four  amplitude-weighted  point  sources,  the 
unnormalized  directivity  function  is: 

D(f,fx,fy)  =  2Q  cos[7Td(f'x  +  fy)]  -  2Q  cos(7rd(f'x  -  fy)]  (11) 

where  =  u  X  =  sinGcosy'X  »  f  =  v/X  =  sinGsiny  X  and  the 

direction  cosine  v  =  sin©  sinNk.  This  yields 

D(f,Q,\|/)  =  -4Q  sin|(md'X)  sinGcosy]  sin((7td  X)  sinGsiny]  (12) 


Figure  2.6  Tesseral  Quadrupole  Geometry 
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which  is  the  far  field,  unnormalized  directivity  function 
for  this  quadrupole  as  a  function  of  ©  and  V.  [Ref.  8]  The 
above  development  is  particularly  important  because  this 
type  array  was  investigated  in  the  course  of  this  work. 

C.  GENERAL  QUADRUPOLE 

Park,  in  his  thesis  (Ref.  1),  developed  an  equation  for 
calculation  of  a  four  element  array  Fig.  (2.7).  The 
equation  was  developed  by  adding  complex  pressures  from  four 
sources  to  yield 

p  =  _!_  ej(wt  -  kR  -  kDjSin0sin<I>)  +  ej(wt  -  kR  +  kD2sin0sin<J>  -  (p.,)  n  3> 
R  R  * 

+  rc  j(wt  -  kR  +  kD,sin0cos<J>  -  <p,) 

R 

+  ej(wt  ‘  kR-  *  kD4sin0cos<t>  -  (p.) 

R 

where  [(1  +  A  +  B  +  C)  R]  exp[j(tot  -  kR)]  is  the  maximum  possible 
pressure  amplitude  and  the  distances  D-^  D2  f  D3  ^  D ^  and 
measured  from  the  origins  to  the  particular  element.  The 
directivity  is: 

=  - {[cos(kD,sin0sin<!>)  +  Acos(kD-,sinOsin<T>  -  9,)  (14) 

(1  +  A+B  +  C)  1  * 

+  Bcos(kD3sin0cos<J>  -  <p3)  +  Ccos(kD4sin0cos<I>  +  <p4)]“ 

+  (sin(kD|SinOsin<I>)  -  Asin(kD,sin0sin<I>  -  9,) 

-  Bsin(kD^sinOcos<I>  -  9,)  +  Csin(kI)4sin0cos(P  +  <p4>]~ }  '2 
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Figure  2.7  Quadrupole  Geometry 

If  D}  =  D2  =  D3  =  D4  then  comparison  of  Figs.  (2.7)  and 
(2.6)  shows  that  the  configuration  is  simply  the  Tesseral 
quadrupole  of  Fig.  (2.6)  where  ^ has  been  shifted  45°. 

Additionally,  the  axial  quadrupole  of  Fig.  (2.5)  is 
obtained  using  Park's  program  by  placing  source  B  at  an 
infinitesimally  small  distance  from  the  array  center 
(i.e.,  D3  =  0.00001  meters  for  any  A  )  with  amplitude 

weighting  equal  to  twice  sources  1  and  A  and  source  C 
amplitude  weighted  to  zero,  which  removes  its  contribution 


2  4 


from  the  array.  This  axial  quadrupole  can  also  be  obtained 
by  placing  sources  B  and  C  at  the  array  center  and  amplitude 
weighting  equally. 

These  equations  were  utilized  to  provide  verification 
data  for  the  laboratory  experiments.  The  three  element 
array  data  were  verified  by  amplitude  weighting  element  C  of 
Eq.  (14)  to  zero.  Dipole  data  were  verified  by  amplitude 
weighting  elements  B  and  C  of  Eq.  (14)  to  zero. 
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III.  EXPERIMENTAL  PROCEDURE  AND  DEVELOPMENT 


A.  BEAM  PATTERN  MEASURING  SYSTEM  OVERVIEW 

The  beam  pattern  measuring  system  is  an  instrumentation 
package  which  was  designed  to  power  and  rotate  an  acoustic 
array  and  then  receive  the  signal,  compute  the  pressure 
distribution  of  the  array  and  display  this  data  in  both  a 
tabulated  format  and  plot  of  sound  pressure  level  in 
decibels  versus  degrees  of  rotation.  All  data  were  taken  in 
an  anechoic  chamber  to  minimize  surface  interference. 

The  system  used  for  measuring  beam  patterns  of  different 
arrays  is  divided  into  two  major  groups.  These  are  the 
signal  generation  and  input  section  and  the  automated 
advancement,  data  collection  and  display  section.  An 
overview  of  these  is  now  presented  along  with  a  block 
diagram  of  the  system.  Fig.  (3.1).  A  detailed  description 
of  each  piece  of  equipment  is  contained  in  Appendix  A. 

The  Hewlett-Packard  (HP)  Function  Generator  Model  3314A 
was  utilized  for  all  signal  generation.  The  signal  passed 
through  an  individual  HP-467A  Power  Amplifier  to  each  of  the 
University  Sound  Model  ID-30  acoustic  driver  units.  Setting 
the  output  levels  of  the  driver  units  was  accomplished 
utilizing  the  variable  gain  control  on  the  power  amplifiers. 
This  constituted  the  signal  generation  and  input  section. 
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Figure  3 . 1  System  Diagram 
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The  receiving  and  processing  section  was  controlled  by 

the  HP-86B  computer.  It  directly  controlled  the  array 

rotation  and  the  reading  of  all  measurements  of  peripheral 

instruments  via  the  Hewlett-Packard  Interface  Bus  (HPIB) . 

The  computer  sent  clockwise  stepping  commands  to  the  motor- 

power/direction  controller  via  the  HP-3421A  Data 

Acquisition/Control  (DAC)  to  rotate  the  array  in  stepped 

increments  of  about  one  degree.  Following  each  advancement 

of  the  array,  the  computer  read  both  the  resistance  output 

of  the  array  rotation  measurement  potentiometer  and  the 

voltage  output  of  the  General  Radio  Model  1962  Microphone. 

The  array  rotation  potentiometer  was  a  0  to  10  Rohm 

potentiometer  mounted  on  top  of  the  rotation  motor  and 

measured  the  angle  of  rotation  of  the  array  (0  to  10  Kohms 

corresponded  to  0°  to  360°).  These  resistance  and  voltage 

readings  were  then  stored  in  the  computer  for  use  in  the 

data  analysis  and  plotting  routines.  Once  the  data  was 

o 

taken  for  360  of  rotation,  it  was  used  to  compute  the 
rotation  angle  and  its  corresponding  decibel  (dB)  level  for 
each  datum.  These  values  were  then  output  to  the  computer 

monitor  and  plotted  in  the  form  of  a  polar  plot  of  dB  versus 

degrees  of  rotation. 

B.  BEAM  PATTERN  PROGRAM 

The  Beam  Pattern  program  listed  in  Appendix  B,  consists 
of  a  series  of  sub  routines.  The  program  starts  by 
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providing  addresses  for  the  two  multimeters  and  the  data 
acquisition/control  unit.  Next  the  initialization 
subroutine  resets  and  initializes  these  three  instruments. 

Once  the  equipment  is  ready  to  operate,  the  incrementing 
subroutine  is  called.  This  is  an  iterative  and  interactive 
routine  used  to  rotate  the  array  and  measure  the  voltage  of 
the  sensed  pressure  field  and  the  corresponding  resistance 
reading  which  is  proportional  to  the  angular  position  of  the 
array.  Up  to  360  data  points  can  be  collected  per  data  run. 
After  the  initial  datum  is  read,  the  array  rotator  motor  is 
energized  for  0.5  milliseconds  to  increment  the  array 
azimuth.  The  array  is  then  allowed  to  coast  down  for  4 
seconds  prior  to  reading  the  voltage  and  resistance.  This 
selection  of  rotate  and  coast  down  times  is  discussed  later 
in  the  chapter.  Once  the  resistance  and  voltage  readings 
are  taken,  they  are  stored  in  an  array  for  later  transfer  to 
a  data  disk  and  for  use  in  the  plot  subroutine. 

Finally  the  plot  subroutine  calculates  the  relative 
decibel  level  of  the  beam  pattern  based  on  the  measured 
voltages.  It  then  takes  the  resistance  reading  and  converts 
it  to  the  angle  of  rotation.  The  calculated  decibel  levels 
are  then  normalized  to  zero  decibels  based  on  the  highest 
received  voltage  reading.  Finally,  the  data  points  are 
plotted  on  a  polar  plot  of  relative  sound  pressure  level  in 
decibel  versus  angle  of  rotation. 
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All  raw  data  is  stored  on  a  data  disk  after  a  complete 
data  run.  Each  data  run  is  stored  by  run  number,  as  shown 
in  Table  II.  It  can  be  retrieved  and  printed  at  the 
printer  as  received  voltage  and  resistance  for  each  datum 
point  or  printed  and  plotted  as  decibel  versus  angle  of 
rotation . 

The  physical  weight  of  an  array  was  a  major  variable  in 
this  program  with  regard  to  the  selection  of  the  number  of 
data  points  needed.  Arrays  weighing  less  than  4  pounds 
required  360  data  points  for  a  full  360°  of  rotation.  The 
number  of  data  points  required  for  larger  arrays  was 
decreased  as  array  weight  increased.  For  the  four  element 
array,  weight  was  substantial  enough  (19  pounds)  to  require 
only  317  data  points  for  a  full  360°  rotation.  This  was 
caused  by  the  fact  that  the  heavier  arrays  would  coast-down 
for  a  larger  distance  after  the  0.5  millisecond  motor-on 
period . 

Another  contributor  to  variation  in  the  required  number 
of  data  points  was  array  balance.  When  utilizing  a  three 
element  array,  the  array  was  sufficiently  out  of  balance  to 
require  a  greater  number  of  data  points  for  full  rotation. 
This  was  attributed  to  minor  binding  of  the  rotator  shaft  as 
the  array  coasted  down. 

Selection  of  motor  run  time  and  coast-down  time  was 
based  on  experimental  investigation.  A  0.5  millisecond 
motor  on  time  provided  1 °  rotation  regardless  of  coast-down 
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time  with  a  light  weight  array.  The  4  second  coast-down  was 
selected  because  smaller  periods  did  not  allow  the  array  to 
stop  and  the  resistance  multimeter  to  stabilize  prior  to  it 
being  read.  Motion  at  the  variable  potentiometer  when  the 
resistance  values  were  read  yielded  erratic  reading  on  the 
multimeter  due  to  electromagnetic  interference  generated  by 
the  motor. 

C.  EXPERIMENTAL  APPARATUS  DEVELOPMENT 

The  experimental  apparatus  was  designed  to  optimize  a 
thorough  search  of  two,  three  and  four  element  arrays  and  to 
establish  their  respective  beam  patterns.  Previously, 
several  experimental  investigations  have  been  conducted 
utilizing  a  two  microphone  array  with  a  single  source.  It 
was,  therefore,  decided  that  reciprocity  would  be  invoked 
and  the  array  would  consist  of  multiple  ID-30  driver  units 
(Appendix  A)  with  a  single  microphone  utilized  for  sensing 
the  pressure  field. 

Six  ID-30  driver  units  were  checked  individually  in  the 
anechoic  chamber  for  their  beam  pattern  impedence, 
resistance  and  reactance  characteristics.  This  data  was 
utilized  to  select  the  best  matched  set  of  four  driver 
units.  All  drivers  were  examined  between  5  Hz  and  3000  Hz 
in  an  effort  to  locate  and  identify  the  resonances  in  the 
individual  speakers.  Fig.  (3.2)  shows  a  representative 
resistance  curve  for  these  drivers  between  the  above 
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RESISTANCE  VS 
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Figure  3.2  Typical  ID-30  Resistance  Plot 


frequencies.  It  was  desired  to  conduct  all  analysis  at 
frequencies  in  relatively  flat  areas  away  from  resonances. 
This  was  done  to  assure  phase  stability  since  reactance  is 
not  varying.  The  speakers  utilized  in  this  investigation 
had  resonances  as  per  Table  I. 

TABLE  I 

ID-30  DRIVER  UNIT  RESONANT  FREQUENCIES 
Speaker  Number  :  First  Resonance  (Hz) _ :  Second  Resonance  (Hz) 


1 

689 

2290 

2 

660 

2225 

3 

780 

2150 

4 

798 

2129 

5 

740 

2253 

6 

778 

2125 

Matched  set  number  1  consisted  of  speakers  1  and  2  and  set 
number  2  of  speakers  3  and  4.  All  array  beam  patterns  were 
run  at  300  Hz  and  960  Hz  which  fell  in  the  flat  areas  of  the 
impedance  curves. 

The  array  existed  in  two  individual  configurations,  a 
linear  and  a  crossed  axis  system  with  variable  tilt  (Fig. 
(3.3)).  The  linear  system  was  used  for  investigation  of 
doublets,  some  triplets  and  quadruplets  (in  the  form  of  an 
axial  quadrupole) .  The  crossed  axis  array  proved  to  be 
extremely  versatile  and  useful,  providing  four  arms,  each 
70  centimeters  in  length,  for  holding  the  drivers. 


Additionally,  it  provided  a  means  to  vary  the  angle©  a  full 
90°.  Each  driver  was  equipped  with  a  harness  assembly  which 
mounted  to  the  array  rod  in  such  a  manner  that  the  driver 
face  was  located  2.5  centimeters  off  the  axis.  This  did  not 
prove  to  be  a  matter  of  concern  since  all  pressure 
measurements  were  made  in  the  far  field. 


Figure  3.3  Array  Configurations 
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A  scattering  and  shading  problem  of  small  magnitude  was 
encountered  because  of  the  large  center  piece  on  the  crossed 
axis  array.  This  scattering  problem  was  largely  resolved  by 
wrapping  the  center  piece  with  fiberglass. 

The  General  Radio  Model  1962  microphone  was  located  3.2 
meters  from  the  center  of  the  array  placing  it  in  the  far 
field.  It  was  mounted  at  the  same  height  in  the  anechoic 
chamber  as  the  array  center.  Changes  in  tilt  angle  theta 
were  accomplished  by  tilting  the  array. 

The  outputs  of  the  ID-30  driver  units  were  balanced 
prior  to  each  data  run.  It  was  done  by  aligning  the  array 
in  the  vertical  plane  with  the  drivers  facing  the 
microphone.  Then  the  output  of  each  driver  was  adjusted 
utilizing  the  variable  gain  control  on  the  power  amplifier. 
This  allowed  all  the  drivers  to  be  driven  within  a  received 
voltage  of  0.002  VAC  of  each  other  out  of  a  typical  level  of 
1  Volt. 

Output  signals  from  the  HP-467A  Power  Amplifiers  were 
monitored  on  an  oscilloscope  as  was  the  received  signal  from 
the  microphone.  Because  the  signals  were  clean  and  without 
any  substantial  noise,  filters  were  not  utilized  in  the 
microphone  circuit. 
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IV.  DATA  ANALYSIS  AND  VERIFICATION 


After  establishing  the  ability  to  accurately  rotate  an 
acoustic  array  while  collecting  and  displaying  sensed 
pressure  levels  under  computer  control,  it  was  necessary  to 
ensure  the  data  were  accurate.  This  could  be  done  by  taking 
a  set  of  parameters  and  obtaining  both  theoretical  and 
experimental  beam  patterns  for  them.  Verification  was 
conducted  utilizing  the  results  of  a  FORTRAN  program 
developed  by  Park  (Ref.  1)  which  yields  a  variety  of  useful 
plots  with  the  help  of  DISSPLA. 

A.  THEORETICAL  BEAM  PATTERN  PREDICTION 

With  the  program  developed  by  Park,  quick  beam  pattern 
plots  could  be  made  to  compare  with  experimental  patterns. 
As  an  example.  Fig.  (4.1)  depicts  a  three  dimensional 
representation  of  a  dipole.  Clearly,  this  particular  is  not 
of  much  use. 

The  three  dimensional  plot  was  sectored  for  greater 
utility.  Fig.  (4.2)  shows  a  sector  representation  for  a 
four  element  array  with  kD^  =  kD2  =  3.52  and  kD^  =  kD^  =  0. 
Note  that  the  kD^  utilized  in  the  remainder  of  the  thesis 
related  directly  to  the  spacing  from  the  center  of  the  array 
to  the  element  as  related  to  Fig.  (2.7).  Refer  to  Fig. 
(2.7).  The  Z-axis  of  the  plot  is  beam  pattern  (dB) ,  the  Y- 
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DIRECTIVITY  PATTERN 


kD1=kD2=  3.52 
kD3=kD4=  0 

Dl=02  D2=0 2 
D3=0.  D4=0. 

VIEW  ANGLE=  20,  20 
P2=10O  P3=0  P4=  0 
A=l.  B=0.  C=0. 


Figure  4.1  Three  Dimensional  Dipole  Representation 
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Figure  4. 


kD1=kD2=kD3=kD4=  10.55 

£>1=0.6  D2=  0.6 
£3=0.6  D4=0 .6 
VIEW  ANGLE=+20,+60 
P2=  0  P3=  0  P4=  o 
A=l.  B=l.  C=l. 


2  Sectored  Three  Dimensional  Plot  (  4  Element  ) 
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axis  psi  (degrees),  and  the  X-axis  theta  (degrees). 
Directivity  is  scaled  from  20  to  -40  dB  (with  0  dB  being  the 
highest  achievable  value  because  of  normalization) ,  psi  from 
180°  to  360°,  and  theta  from  0°  to  -90°.  This  plot  appears 
to  be  very  busy  but  proves  to  be  quite  useful.  It  is  easy 
to  evaluate  this  sectioned  representation  by  following  the 
contour  line  of  the  directivity  above  the  line  theta  equals 
0°  for  psi  from  180°  to  360°.  This  provides  the  same  dB 
levels  at  180°  <  psi  <  360°  as  the  polar  plot  representation 
plotted  in  Fig.  (4.3). 

B.  EXPERIMENTAL  DATA-SINGLE  ELEMENT  DIRECTIONALITY 

As  stated  in  Chapter  III,  testing  was  conducted  at  300 
Hz  and  960  Hz  which  were  in  the  flat  areas  of  the  drivers' 
resistance  curves,  away  from  their  resonance  peaks.  This 
was  done  to  assure  phase  stability  between  drivers  since  the 
reactance  is  not  varying  (large  phase  shifts  do  not  occur 
except  near  resonance) . 

Initial  data  runs  were  conducted  with  a  single  ID-30 
acoustic  driver  to  establish  the  directionality  of  each  of 
the  drivers.  The  drivers  were  rotated  360°  at  three 
separate  tilt  angles  of  the  driver  face,  ©  =  0°,  45°  and  90. 
A  representative  graph  of  the  directionality  is  contained 
in  Fig.  (4.4)  which  was  made  with  the  driver  at  a  tilt  angle 
of  0°  and  pointing  at  the  microphone  prior  to  rotation.  The 
dotted  line  is  the  laboratory  beam  pattern.  Note  that  there 
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PHI  VS.  DIRECTIVITY 


0 


THETA=  0  DEORT^ 

kD1=kD2=kr)3=kD4- 


Ul=0.6  D2=0.6 
D3=0.6  D4=0  6 
P2=  0  P3=  0  P4= 
A=l.  B=l.  C=l. 


10.55 

0 


Figure  4.3  Polar  Plot  of  a  Four  Element  Array 
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Figure  4.4  ID-30  Acoustic  Driver  Directionality 


is  a  1  to  2  dB  drop  below  the  0  dB  level  as  the  driver 
points  away  from  the  receiver  which  shows  the  driver  is  not 
completely  omnidirectional.  It  should  also  be  noted  that 
the  program's  normalization  routine  sets  the  highest 
received  microphone  voltage  level  to  0  dB  and  normalizes  all 
other  voltage  levels  in  the  beam  pattern  to  this  level. 
This  was  done  to  provide  greater  ease  in  plotting  of  beam 
patterns  so  as  to  produce  a  uniformity  between  plots. 

C.  EXPERIMENTAL  AND  THEORETICAL  DATA  -  VARIOUS  ARRAY 

CONFIGURATIONS 

An  investigation  of  two,  three  and  four  element  arrays 
was  conducted.  For  each  of  these  configurations,  several 
different  combinations  of  spacing  and  amplitude  and  phase 
shading  were  looked  at.  Table  2  at  the  end  of  this  chapter, 
is  a  list  of  these  variations  for  which  beam  patterns  were 
obtained.  Note  that  phase  shading  was  limited  to  either  $ 
or  180°.  In  viewing  the  figures,  the  term  "x-pole"  is  used 
loosely  to  represent  the  term  "x-let"  (i.e.,  quadrupole  for 
quadruplet) . 

For  the  two  element  array,  beam  patterns  were  run  in  the 
laboratory  for  both- in- phase  and  180°  out-of-phase  doublets. 
An  example  of  the  out-of-phase  doublet  beam  pattern  obtained 
in  the  anechoic  chamber  is  shown  in  Fig.  (4.5).  (Note  that 
in  this  figure  kd  is  the  quantity  2kDf  when  Dj  =D2*) 

Its  theoretical  verification,  obtained  from  Park's  program 
is  shown  in  Fig.  (4.6).  A  contributing  factor  to  their 
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THETA  VS.  DIRECTIVITY 

0 


PHI-  90  DEGREE 

kD1=kD2=  5.276  ,  kD3=kD4=  0 

Dl=D2=  0.5 
D3=D4=  0 

P2=180  P3-  0  P4=  0 
A=l.  D---0.  C=0. 


Figure  4.6  Theoretical  Dipole  Plot 
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small  differences,  is  the  slight  directionality  associated 
with  the  acoustic  driver.  The  slight  offset  of  the 
laboratory  plot  nulls  is  caused  by  the  array  not  being 
exactly  aligned  to  broadside  at  0  Kohm  on  the  potentiometer. 
The  90°  rotation  between  the  two  plots  is  caused  when  theta 
is  fixed  at  90°  vice  0°  for  theoretical  plots.  Because  of 
the  narrowness  of  the  nulls,  laboratory  data  would  require  a 
sample  exactly  with  the  array  broadside  to  achieve  the  null. 
Other  than  the  difference  in  nulls,  there  is  exceptionally 
good  agreement  between  experimental  and  theoretical  beam 
patterns . 

In  the  investigation  of  three  element  arrays,  the  arrays 
were  configured  as  equilateral  triangles,  isosceles 
triangles  and  as  a  linear  array.  For  phase  shading,  three 
cases  were  looked  at:  all  elements  in  phase,  the  center 
element  180  out-  of-phase,  and  one  end  element  placed  180 
out-of- phase .  These  arrays  were  also  subjected  to  various 
amplitude  shading  (the  axial  quadruplet  composed  of  three 
elements  with  the  center  element  amplitude  shaded  to  twice 
the  value  of  each  of  the  end  elements  will  be  discussed  in 
the  quadruplet  section) .  Experimental  data  for  a  three 
element  isosceles  triangle  with  all  elements  in  phase  and 
kD^  =  kD2  =  2kD3  =  5.28  are  shown  in  Fig.  (4.7)  with  its 
theoretical  verification  shown  in  Fig.  (4.8).  The  deep 
theoretical  nulls  seen  on  Fig.  (4.8)  are  not  present  in  the 
laboratory  data.  This  can  be  attributed  to  a  small  variance 


45 
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THETA  VS.  DIRECTIVITY 


PHI=  90  DEGREE 

kD1=kD2=2kD3=  5.276 

Dl-03  D2=0.3 
D3=0.I5  D4=0. 

P2=  0  P3=  0  P4=  0 
A-l.  B=l.  00. 


Figure  4.8  Three  Element  Isosceles  Plot 
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in  the  equal  amplitude  shading,  a  small  amount  of  scatter 
from  the  large  array  center  joint,  and  a  slight 
directionality  of  the  individual  acoustic  drivers. 

The  beam  patterns  of  a  triplet  in  the  form  of  an 
isosceles  triangle  are  shown  in  Figs.  (4.9),  (4.10),  and 

(4.11)  to  show  the  effects  of  tilting  the  array  out  of  the 

horizontal  to  angles  of  theta  equal  to  30°  ,  60°  and  9(P 

respectively.  The  plane  of  the  elements  was  initially 

perpendicular  to  the  x-y  plane,  making  theta,  the  angle 

o 

between  the  z-axis  and  the  array  plane  equal  to  0  .  The 

array  plane  was  then  tilted  30°  about  the  x-axis  and  rotated 
360°  about  the  z-axis  to  obtain  the  data  for  Fig.  (4.9). 
The  array  plane  was  then  tilted  in  30  increments  to  obtain 

Figs.  (4.10)  and  (4.11).  For  this  triplet  case,  kD  =  kD  2  = 

2kD3  =  5.28  with  kD3  sensing  the  elevation  change  (Fig. 

(4.12) )  . 
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Figure  4.11  Beam  Pattern  for  90  Tilt  in  Theta 


Figure  4.12  Array  Configuration  for  Tilt  Theta 

These  variations  in  tilt  angle  (theta)  show  that  the  3^ 
tilt  has  a  diminished  dB  level  for  all  lobes  except  the 
broad  lobe  coming  into  and  exiting  broadside.  This  is  due, 
in  large  part,  to  the  directionality  of  the  drivers,  because 
they  are  facing  the  microphone  to  a  greater  extent  at  31? 
than  at  90°  of  tilt  and  mutual  shading  exists. 
Additionally,  for  these  various  tilt  angles,  the  Z-axis  is 
rotating  about  the  local  vertical,  as  explained  in  section  D 
of  this  chapter. 

The  last  major  configuration  investigated  was  the 
quadruplet.  It  was  looked  at  both  as  an  axial  quadruplet. 
Fig.  (2.5),  and  as  a  planar  array,  Fig.  (2.7).  An 
experimental  beam  pattern  for  a  planar  array  quadruplet  with 

^D1  =  kD2  =  kD3  =  kD4  =  5.26  and  all  elements  in  phase  is 
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shown  in  Fig.  (4.13)  along  with  a  theoretical  beam  pattern, 
Fig.  (4.14).  There  is  good  agreement  between  the  two  plots. 

A  plot  of  an  axial  quadruplet  represented  by  an 
amplitude  shaded  triplet,  obtained  in  the  laboratory  is 
shown  in  Fig.  (4.15).  Comparing  Fig.  (4.15)  to  Fig.  (4.13), 
it  is  seen  that  when  two  elements  of  the  planar  array  are  in 
line,  the  arrival  from  the  more  distant  source  is  weakened 
by  the  shadowing  effect  of  the  nearer  source.  Therefore,  the 
near  source  can  be  emphasised  by  as  much  as  one  dB  when 
compared  to  the  far  source.  Additional  plots  of  various 
configurations  are  contained  in  Appendix  C,  where  good 
agreement  continues  to  be  seen. 

D.  EFFECTS  OF  ARRAY  TILT,  THETA 

Referring  back  to  Fig.  (4.2),  it  is  again  noted  that 
contour  lines  are  directly  related  to  polar  plots  when 
either  theta  or  psi  is  fixed.  However,  the  countour 
representation  is  not  as  useful  when  theta  and  psi  are  both 
changing  with  rotation  of  the  array  as  in  the  case  where  the 
array  has  been  tilted  in  theta.  In  this  case  the  polar 
curve  is  found  by  superimposing  a  trajectory  which  is  the 
specific  locus  of  @  and  ^  in  the  ©  -'ty  plane  onto  the 
contoured  surface. 

For  this  work,  the  microphone  was  fixed  in  the  far  field 
at  a  distance  of  3.2  meters  and  level  with  the  center  joint 
of  the  array.  Tilting  the  array  in  angle  theta  caused  an 
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theta  vs.  directivity 


0 


PHI=  90  DEGREE 

kDl=kD2=kD3=kD4=  5-275 

Dl=05  D2-0  5 
D3=0.5  D4=0  5 
P2=  0  P3=  0  P4=  0 
A=l.  B-0.  C=0. 


Figure  4.14  Theoretical  Planar  Quadruplet  Plot 
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R-3.20  M  kD1=kD2-kD3=  4.713 

PHASE-O.  100.  0  Deg  AMPL-  0.  491  Vac  SHADE-1.2.  1 


offset  between  the  Z-axis  and  the  local  vertical  about  which 


the  array  rotated.  This  can  be  seen  in  Fig.  (4.16).  In  the 
left  sketch,  the  array  elements  are  aligned  in  the  standard 
rectangular  coordinate  system,  with  the  elements  in  the  x-y 
plane  and  rotating  about  the  z-axis.  However,  when  the 
plane  of  the  array  is  tilted  30°  about  the  x-axis,  Y'  =  Y+3(P 
and  z  '  =  z+30° .  A  90°  rotation  about  the  z-axis  to  21$ 
yields  z  '  =  90°  -  30°  =  60°  .  This  means  that  a  36l^ 
rotation  about  the  z-axis  causes  Yl  and  Zl  to  vary  from  12$ 
to  60°  and  back  to  120°  .  This  is  due  to  the  fact  that  the 
microphone  is  fixed  in  space  in  the  far-field  rather  than 
tilting  with  the  plane  of  the  array. 


Figure  4.16  Z-axis  Offset  Effects 

This  problem  does  not  occur  in  the  classical  representation 
where  the  microphone  is  fixed  with  respect  to  the  axis  of 
rotation  of  the  array. 

E.  ERRORS 

The  sources  of  errors  in  the  beam  pattern  measuring 
system  were  examined  both  mathematically  and  empirically. 
Sources  of  error  for  direct  error  analysis  were  based  on 
multimeter  values  and  power  amplifier  settings.  With  output 
reading  from  the  multimeters  and  power  amplifier  exhibiting 
precision  to  the  third  decimal  place,  instrumentation  errors 
from  these  sources  are  less  than  1%  and  are  not  major 
contributing  errors  to  measured  beam  patterns. 


58 


The  major  source  of  error  in  the  system  was  determined 

empirically  to  be  the  acoustic  driver  positioning  on  the 

array.  The  effects  of  a  small  change  in  kd  are  shown  in 

Figs.  (4.17),  (4.18),  and  (4.19).  Here,  it  is  seen  that 

small  changes  in  kD^  prove  to  have  a  large  effect  on  the 

size  of  the  small  lobe  at  ±  90°  and  0°  and  180°  ,  thus 

yielding  the  largest  error.  A  change  of  ±  5%  in  kd  effected 

the  structure  by  up  to  6  dB  as  seen  from  the  behavior  of  the 

o 

small  lobe  at  ±_  90  .  This  change  relates  to  a  positioning 
error  of  about  ±  3  centimeters  and  presented  a  small  problem 
with  respect  to  exact  positioning  of  the  speakers. 

F.  VALIDATION  OF  RESULTS 

The  laboratory  beam  patterns  prove  to  be  in  good 
agreement  with  theoretical  beam  patterns.  Small  uniform 
beam  pattern  amplitude  differences  are  merely  a  consequence 
of  the  difference  in  normalization  methods  employed  in  the 
theoretical  and  experimental  beam  pattern  processing.  The 
theoretical  program  normalized  the  directivity  to  the  sum  of 
the  amplitude  weights,  whereas  the  laboratory  program 
normalized  the  beam  pattern  to  the  highest  received  voltage. 

An  independent  verification  of  the  experimental  program 
was  performed  by  a  thesis  group  under  the  auspices  of 
Professor  Don  Walters  of  the  Naval  Postgraduate  School 
Physics  Department,  who  asked  to  use  this  experimental 
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THETA  VS.  DIRECTIVITY 


0 


PH1=  60  DEGREE 

kD1=kD2=kD3=kD4=  4.75 

Dl=0.5  D2=05 
D3=05  D4=0.5 
P2=  0  P3=  0  P4=  0 
A=l.  B=l.  C=l. 


Figure  4.17  kDi=  4.75  for  Error  Computation 
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THETA  VS.  DIRECTIVITY 

0 


PHI=  60  DEGREE 

kD1=kD2=kD3=kD4=  5.0 

U1=0  5  D2=0.5 
D 3=0.5  D4=0.5 
P2=  0  P3=  0  P4=  0 
A=l.  B=l.  C=l. 


Figure  4.18  kD^-  5.0  for  Error  Computation 
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THETA  VS.  DIRECTIVITY 

0 


PHI=  60  DEGREE 

kDl=kD2=kD3=kD 4=  5.25 

Dl=0.5  D2=05 
D3=05  D4=0.5 
P2=  0  P3=  0  P4=  0 
A=l.  B=l.  C=l. 


Figure  4.19  kD-[=  5.25  for  Error  Cortputation 
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The  laboratory  system  was  placed 


system 
in  use 
for  a  5 
further 


for  related  work. 

by  this  research  group  in  the  selection  of  drivers 
x  5  array  currently  under  investigation  and  will  be 
utilized  to  establish  beam  patterns  for  this  array. 


TABLE  II 


ARRAY  CONFIGURATIONS 


Run  No/ 

Amplitude  Spacing 

No. 

Preq 

Phase  Tilt 

Shad ing 

from  Center 

Elements  ; 

:  (Hz) 

t 

(deg)  i  (deg)  i 

(VAC) 

i  (CN) 

1 

300 

0 

1 

300 

45 

1 

300 

90 

1 

960 

0 

1 

960 

45 

1 

960 

90 

2 

300 

0,0 

0 

30 

2 

300 

0,180 

0 

30 

2 

960 

0,0 

0 

30 

2 

960 

0,180 

0 

30 

2 

960 

0,0 

90 

30 

2 

960 

0,180 

90 

30 

2 

300 

0,0 

90 

30 

2 

300 

0,180 

90 

30 

2 

300 

0,0 

45 

30 

2 

300 

0,180 

45 

30 

2 

960 

0,0 

45 

30 

2 

960 

0,180 

45 

30 

4 

960 

0,0, 0,0 

0 

30 

4 

300 

0,0, 0,0 

0 

30 

4 

960 

0 

,180,0,180 

0 

30 

4 

300 

0 

,180,0,180 

0 

30 

4 

960 

0,0, 0,0 

30 

30 

4 

300 

0, 0,0, 0 

30 

30 

Quadll/4 

300 

0 

,180,3,180 

30 

30 

Quadl 2/ 4 

300 

0 

,180,0,180 

60 

30 

Quadl3/4 

960 

0 

,180,0,180 

60 

30 

Quadl 4/ 4 

960 

0,0, 0,0 

60 

30 

Quadl5/4 

300 

0, 3,0,3 

60 

30 

Quadl 6/4 

960 

0,0, 0,0 

90 

30 

Quadl7/4 

960 

0 

, 180,0,180 

90 

30 

Quadl8/4 

960 

0 

,180,0,0 

0 

30,30,30,0 

Tr i 1/3 

960 

0,0,0 

0 

30 

Tri2/3 

300 

0,0,0 

0 

30 

Tri 3/3 

300 

0,180,0 

0 

30 

Tr i 4/3 

960 

0,180,0 

0 

30 

Tr i 5/3 

960 

0,0,180 

0 

30 

Tr i6/3 

960 

0,3,180 

30 

30 

Tri^/3 

960 

0,180,0 

30 

30 

Tri8/3 

960 

3,3,3 

30 

30 

Tri9/3 

960 

0,3,3 

30 

30,15,30 

Tr i 10/3 

960 

0,3,3 

0 

30.15.30 

Trill/3 

960 

0,180,0 

0 

30,15,30 

Tr il2/3 

960 

0,180,0 

30 

30,15,30 

Tr 1131/3 

960 

0,180,0 

60 

30,15,30 

Tr i 14/3 

960 

3,0,3 

60 

30,15,30 

Tr i 15/3 

300 

0,0,0 

60 

30,15,30 

Tr i 16/3 

960 

0,0,0 

90 

30,15,30 

Tr 117/3 

960 

0,180,0 

90 

30,15,30 

Tr i 18/3 

960 

0,0,180 

90 

30,15,30 

Tri20/3 

960 

0,0,0 

0 

1,1,1 

30,0,30 

Tr i 2 1 / 3 

960 

0,180,3 

0 

1,1,1 

30,3,30 

Tr i22/3 

960 

0,0,3 

0 

1,2,1 

30,0,30 

Tr i23/3 

960 

0,180,3 

0 

1,2,1 

30,0,30 

Tr i2  4/3 

960 

0,0,180 

0 

1,2,1 

30,0,30 

Tr i25/3 

960 

0,3,3 

0 

1,2,1 

30,30,30 

Tr  i 2  6/3 

960 

0,0,0 

30 

1,2,1 

30,30,30 

Tr i 27/3 

960 

0,0,3 

60 

1,2,1 

30,30,30 

Tr i28/3 

960 

0,0,0 

90 

1,2,1 

30,30,30 

Tr i 30/3 

960 

0,0,0 

90 

1,4,1 

30,30,30 

Tr i 31/3 

960 

0,130,0 

90 

1,4,1 

30,30,30 

Tr i 32/3 

960 

0,0,3 

0 

1,2,1 

20,20,20 

Tr i 33/3 

960 

0,180,3 

0 

1,2,1 

20,20,20 

Tr i 3  4/3 

960 

0,180 , 0 

0 

1,2,1 

20,10,20 

Tr i36/3 

960 

0,3,3 

0 

1,2,1 

20,10,20 

Tri 37/3 

960 

0,3,0 

0 

1,2,1 

10,10,13 

Tri 38/3 

960 

0,180,0 

0 

1,2,1 

10,10,10 

Tr i 39/3 

960 

0,180,3 

90 

1,1,1 

20,30,20 

Tr i40/3 

963 

0,180,0 

90 

1,2,1 

17.9,17.9,17.9 

Tr i 41/3 

960 

0,0,3 

90 

1,2,1 

17.9,17.9,17.9 

Tr i43/3 

960 

3,0,0 

90 

1,2,1 

26.8,26.8,26.8 

Tr i 4  4/3 

960 

0,180,0 

90 

1,2,1 

26.8,26.8,26.8 

Note:  (1)  Additional  data  runs  were  made  in  some 

of  the  same  configurations  with  fiberglass 
wrapped  around  the  centerpiece  to  identify 
the  effects  of  scattering. 

(2)  Additional  data  runs  were  made  in  some 
of  the  same  configurations  to  verify 
repeatability . 
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V.  CONCLUSIONS  M2  RECOMMENDATIONS 


The  doublet  forms  the  basis  for  countless  sonar  and 
radar  arrays  currently  utilized  in  both  military  and 
civilian  applications.  It  is  the  ability  to  control  the 
amplitudes  and  phases  of  the  element  of  an  array  which 
allows  the  tailoring  of  the  array  radiation  pattern  to 
achieve  the  desired  element  of  interest  and  promotes  the 
maximization  of  array  performance  by  increasing  signal, 
decreasing  noise  or  both. 

The  ability  to  rapidly  measure  the  beam  pattern  of  an 
array  with  complete  confidence  in  the  results  is  the  basis 
behind  this  study.  A  laboratory  system  was  developed  to 
allow  collection  of  a  360  point  beam  pattern  in  a  45  to  50 
minute  period.  Software  is  easily  modified  to  allow  for 
collection  of  more  or  less  than  360  data  points.  One  of  the 
major  advantages  of  the  computer-controlled  method  of 
obtaining  a  beam  pattern  is  the  relatively  small  amount  of 
time  required  to  achieve  an  accurate  beam  pattern  based  on  a 
large  data  base. 

Recommendations  for  improving  the  current  laboratory 
system  are  based  on  reducing  electrical  noise  of  the 
rotation  motor  and  reducing  the  size  and  weight  of  the 
array.  First,  the  noise  of  the  rotation  motor  prohibits 
continuous  rotation  while  reading  the  resistance  on  the 
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potentiometer.  Since  a  CW  signal  does  not  require  broadband 
detection,  filtering  could  be  accomplished  by  use  of  a  lock- 
in,  fast  fourier  transform  or  heterodyne  filter.  Second,  the 
array  weight  should  be  reduced  to  lessen  the  drag  on  the 
rotation  motor  and  to  reduce  the  amount  of  coast-down  time 
required  for  heavier  arrays.  Next,  the  array  center  joint 
should  be  reduced  in  size  to  lessen  any  shadowing  effects 
from  the  large,  metal  center  joint.  (Scattering  effects 
seemed  to  be  small  judging  from  the  lack  of  any  effect  on 
the  polar  plots  when  this  joint  was  wrapped  in  fiberglass.) 
Also,  driver  mounts  need  to  be  designed  so  that  the  center 
of  the  driver  face  is  directly  on  the  axis.  The  2 
centimeter  offset  utilized  for  this  investigation  was 
minimized  by  making  all  readings  in  the  far  field.  Last, 
greater  distance  between  the  driver  and  microphone  would 
also  help  to  lessen  the  1  dB  variation  seen  when  an  array 
was  in  an  endfire  orientation.  This  only  affected  the  depth 
of  the  nulls. 
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APPENDIX  A 


EQUIPMENT  DESCRIPTION 

1.  ANECHOIC  CHAMBER 

All  experimentation  was  conducted  in  the  anechoic 
chamber,  Room  019,  Spanagel  Hall,  Naval  Postgraduate  School. 
The  anechoic  chamber  is  a  99%  echo  free  enclosure  for  sounds 
above  100  Hertz.  Acoustics  research  can  be  conducted  on 
sound  sources,  receivers  and  scatterers  with  minimal 
interference  from  surface  reflections  or  external  noises. 
The  size  of  the  outer  room  in  which  the  anechoic  chamber 
floats  is  37  feet  x  24  feet  x  21  feet  high.  Noise 
transmission  from  outside  is  minimized  by  the  floating  room 
within  a  room  construction.  The  outer  12  inch  concrete- 
walled  room  is  separated  from  the  inner  room  of  concrete- 

block  sides  and  floor  by  a,  2  inch  blanket  of  fiberglass  or 

cork.  Forty  inch  fiberglass  wedges  are  attached  to  all 
surfaces  of  the  inner  room.  These  are  designed  to  trap  and 

absorb  sounds  which  are  incident  upon  them.  A  suspended 

floor  is  constructed  of  a  grid  of  225  wire  cables,  each 
stretched  at  a  tension  of  150  to  200  pounds. 

2.  HP-86B  PERSONAL  COMPUTER 

The  HP-86B  computer,  operating  in  Hewlett-Packard  BASIC 
language,  is  equipped  with  a  64  Kilobyte  (K-byte)  memory 
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with  two  additional  128  K-byte  memory  modules.  The  system 
proved  more  than  sufficient  for  both  program  development  and 
as  an  instrument  controller.  Instrument  control  by  means  of 
the  computer  was  accomplished  with  the  HP-IB  (Hewlett- 
Packard  Interface  Bus)  which  is  a  standard  connection  on 
most  HP  equipment.  This  facilitated  instrument  control  for 
the  various  system  components. 

3.  HP  7470A  GRAPHICS  PLOTTER 

The  HP  7470A  Graphics  Plotter  is  a  vector  plotter  which 
utilizes  HP's  micro-grip  drive  technology,  moving  both  the 
paper  and  the  pen.  Multicolor  plotting  is  available,  but 
there  are  only  two  pen  positions  so  pauses  must  be  used  to 
change  pen  colors.  Only  one  line  of  program  code  is 
required  to  shift  the  output  to  and  from  the  plotter  to  the 
HP-86B  CRT. 

4.  HP-4192A  LOW  FREQUENCY  IMPEDANCE  ANALYZER 

The  HP-4192A  LF  Impedance  Analyzer  is  a  fully  automated 
instrument  designed  to  measure  a  large  range  of  impedance 
related  characteristics.  It  can  be  used  to  measure  eleven 
impedance  parameters:  absolute  value  of  impedance  Z, 
absolute  value  of  admittance  Y,  phase  angle©,  resistance  R, 
reactance  X,  conductance  G,  susceptance  B,  inductance  L, 
capacitance  C,  dissipation  factor  D,  and  quality  factor  Q. 
The  HP-4192A  was  used  to  measure  the  frequence  response. 
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impedance,  resistance  and  reactance  of  the  University  Sound 
ID-30  driver  units  as  discussed  in  section  10  of  this 
appendix . 

5.  HP-3421A  DATA  ACQUISITION/CONTROL  UNIT 

The  HP-3421A  Data  Acquisition/Control  Unit  (DAC)  has  a 
5.5  digit  voltmeter,  a  10  kHz  counter  and  an  HP-IB 
controller  interface.  The  DAC  was  utilized  as  a  computer 
controlled  switching  unit  with  channels  1  and  2  configured 
as  actuators.  This  provided  a  voltage  path  to  control 
relays  in  the  motor  power  and  direction  controller  for 
rotation  of  the  ID-30  speaker  array.  Fig.  (A.l). 

6.  HP-3478A  MULTIMETER 

The  HP-3478A  Multimeter  provides  both  resistance  and 
voltage  measurement  capabilities  with  3.5  to  5.5  digit 
resolution  with  +  0.3  mV  accuracy  on  the  100  Hz  to  20  kHz 
scale.  The  multimeter  was  remotely  controlled  by  the 
computer  via  the  HP-IB.  Voltage  readings  from  the  GR-1962 
microphone  and  resistance  readings  from  the  variable 
potentiometer  on  the  array  rotator  being  read  by  the 
multimeter  and  then  sent  to  the  HP-86B  computer  via  the  HP- 
IB. 


7.  HP-3314A  FUNCTION  GENERATOR 

The  HP-3314A  Function  Generator  is  a  multimode  function 
generator  which  provides  sine,  square,  and  triangular 
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Data  Acquisition  /  Control  Unit 


Figure  A.l.  Data  acquisition  connection  to  motor-power 
controller. 
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functions  as  well  as  any  desired  arbitrary  wave  form  from 
frequencies  of  0.001  Hz  to  19.999  MHz.  It  was  operated  in 
the  free  run  mode  providing  a  continuous  1  Vac  sine  wave  to 
Hp-467A  Amplifier's  inputs. 

8.  HP-465A  AMPLIFIER 

The  HP-465A  Amplifier  is  a  general  purpose  preamplifier 
and  impedance  converter  from  10  megohms  to  50  ohms.  It  has 
a  selectable  gain  of  20  or  40  decibels  and  is  stable  over  a 
frequency  range  of  5  Hz  to  1  MHz.  The  amplifier  was  used  to 
amplify  the  received  signal  from  the  GR-1962  microphone. 

9.  HP-467A  POWER  AMPLIFIER 

The  HP-467A  Power  Amplifier  is  a  10  watt  peak  power 

amplifier  and  a  -20  to  +20  void  dc  power  supply.  It  has  a 

wide  band  width  with  low  dc  drift  from  5  Hz  to  1  MHz.  A 
variable  gain  control  allows  adjustments  of  gain  between 
zero  and  10  with  resolution  better  than  0.1%  of  full  power. 
Four  HP-467A's,  one  for  each  ID-30  driver  unit,  were 
utilized  in  adjusting  power  inputs  from  the  function 
generator . 

10.  UNIVERSITY  SOUND  ID-30  DRIVER  UNITS 

The  University  Sound  ID-30's  are  medium  power  horn 
driver  units  designed  specifically  for  wide  midrange 
response.  Their  tropicalized  voice  coils  are  2  inches  in 
diameter  with  linen  base  molded  phenolic  diaphragms  with  a 
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minimal  16  ohm  impedance.  The  frequency  response  is  85  Hz 
to  7.5  Khz.  They  were  used  as  the  sound  sources. 

11.  GENERAL  RADIO  TYPE  1962  MICROPHONE 

The  GR-1962  microphone  is  a  0.5  inch  diameter  electret 
condenser,  microphone,  providing  a  flat  response  at  zero 
decibels  between  20  Hz  and  15  Khz.  It  experiences  rapid 
drop  off  above  15  Khz  and  is  nonrespons ive  above  21  Khz.  A 
General  Radio  Model  1560-P62  inline  power  supply  was 
connected  to  the  microphone  1560-P41  preamplifier  to 
amplify  the  received  signal. 

12.  HP-6215A  POWER  SUPPLY 

The  HP-6215A  power  supply  provides  zero  to  30  volts  dc 
or  0  to  500  Mamps .  It  was  used  to  provide  switching  voltage 
for  the  Motor  power  relay  in  the  direction  controller. 

13.  MOTOR  POWER  AND  DIRECTION  CONTROLLER 

The  controller  for  the  array  rotator  controlled  power  to 
the  rotator  and  also  the  direction  of  rotation.  With  leads 
to  the  rotator  assembly  potentiometer,  it  also  provided 
resistances  which  corresponded  to  a  given  angle  (0  to  10 
Kohms  translated  into  0  to  360  degrees) .  A  +5  volt  dc 
signal  to  IC-1  via  channel  1  of  the  DAC  provided  on/off 
power  to  the  rotator  and  +5  Volts  to  IC-2  via  channel  zero 
of  the  DAC  provided  clockwise  rotation  of  the  driver 
assembly  unit  (Figure  A. 2).  [Ref.  8] 
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ccw 


a)  Circuit  Diagram  for  Modified  Motor  Power  and  Direction  Controller 


b)  360°  /  1 0  kO  Potentiometer  c)  1 6  Pin  DIP  Relay 


Figure  A. 2.  Circuit  diagram  for  modified  motor-power/  direction 
controller  [Ref.  9.,  p.  57], 


73 


APPENDIX  B 


BEAM  PATTERN  (X)MPUTATION  AND  PLOT  PROGRAM 


10  PRINT  M .  BfcfiM  PATTERN  COMPU  I  AT  I  ON  AND  PLuf  PROGRAM  ***** 

11  FRITH  '* . ****** . ***.*.*♦•♦♦ . * . «" 

1  2  < 

13  i 

23  PRINT  *****  SET_UP  **♦♦*" 

2  1  i 


30 
4  0 
GO 
GO 
6  1 
62 
70 
GO 
33 
100 
1  10 
1  1  1 
1  20 


PRINT  CHRP  <  10  > 

PRINT  "THIS  PROGRAM  CONTROLS  THE  FOLLOWING  APPARATUS  TO  MEASURE,  OTSPLAr,' 
PRINT  " PL 01  mND  SluRE  THE  DIRECTIONAL  CHARACTERISTICS  OF  .A  SINGLET/' 

PRINT  " DOUBLET  ,  TRIPLET,  AND  OUnORUPLET  ARRAY/' 


PRINT  ' 
FRITH  * 
PRINT  ' 
PRINT 
PRINT 

i 

DiSP  ' 


INSTRUMENT  ADDRESS 

HP.  3473  MULTIMETER  (HOLTS)  72  4 
"  HF/ 3473  MULT  IMP  I ER  ‘RES)  723 

"  DnTrt  AQU1S IT  ION 'CONTROL  703 

INITIALIZE  INSTRUMENTS'' 


12  1  i 


130  REM  ***  PRESET  nND  INITIALIZE  INSTRUMENTS*  *  * 

140  0UTPU1  703  ;  " RS "  i  RESET 

1  GO  (HI  RUT  723  ;  T3RST 1 "  •  FREE  RUN  30K0HMS 

ICO  OUTPUT  7-24  ;  "F2RAT  1  "  «  FREE  RUN  AUTOPnNGE  AC  UuLTS 
161  1 

170  LISP  "INITIALIZE  PROGRAM" 

171  i 

ISO  OPTION  BnSE  1 

17)0  COM  nlHLLU  I  G0  )  ,  CALC_ANGLE  (  150),  NAME'S  I  30  ] 

20O  DEG 

2  10  PR  INTER  IS  701 
220  1  NT  EGE  R  1,0, L,N 

2  30  Girl  P'  360  »  /JOLT (  3G0  )  ,GHMS_DATA(  360  )  ,X(  360  )  ,Yv  360  )  ,K(  360  ) 
2  40  ffijLI-0 

2G0  MAGS  STORAGE  IS  ".DATA!" 

260  D I r P  "NUMBER  OF  DATA  POINTS  DESIRED?" 

270  BEEP  10,300 
260  INPUT  H 

2  30  DISP  " WAN  1  NEW  OATA  ( Y / N  ) ? " 

300  PEEP  1O/30O 
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310  INPUT  OS 

320  If  Qi="v"  THEN  GOTO  350 
330  IF  O'f  = " N "  1  MEN  GOTO  700 
340  GOTO  30O 

350  PEN  •••  DEFAULT  PARA.METEFS  *»* 

30 1  i 

323  i  MUST  INSERT  FOR  LINES  390,  4  00,  4  10  AND  810  PRIOR  TO  START  OF  RUN 
303  i 

390  CREATE  "*  INSERT  NAME  OF  RUN  *",360,16 
400  AbSIGM#  I  TO  “*  INSERT  NaME  OF  RUN  *" 

410  NAMES=" •  INSERT  NAME  OF  FIIE  »"  '  FILE  NAME 
420  FULL  SEAL E  =  S00  !  LOCK_IN  MU AG 
4  30  COMU  =  380/ 1 0000  i  C0NUER1  OHMS  TO  DEGREES 
440  5ET..UP:  '  ••*  SET_UP  MEASUREMENT  PARAMETERS  *«* 

4  50  PRINT  CHRS  <  1 0  ) 

4 50  LISP  "MANUALLY  ENTER  DESIRED  FREQUENCY  FOR  SPEAKER  DRIVER" 

500  DISP  "ENTER  0  IF  NO  CHANGE" 

510  DISP  "ENTER  I  TO  CHANGE  TIRE" 

520  EEEP  10,600 
530  INPUT  C  !  CHOICE 
540  IF  C=0  THEN  GOTO  560 
650  IF  C= I  I  HEN  GOTO  750 
560  PRINT  CHPi  (10) 

570  OUTPUT  709  ;"CLS0"  l  ROTATE  CW 

575  PRINT  "  I "  , " uHMS_DAT A  i OHMS  )“,  "VOLT  ( VAC  )" 

590  FOR  1=1  TO  N 

590  OUTPUT  709  ;"CLGI"  i  POWER  TO  MOTOR  ON 
GOO  WAIT  .5  I  MILLISECONDS 

610  OUTPUT  709  ;  "GPN 1  "  I  POWER  TO  MOTOR  OFF 
620  WAIT  4000 

630  ENTEP  723  ;  OHMS_DATA( I  ) 

640  I.N1FR  724  ;  VOLT(I) 

650  PRINT  I  , OHMS  _DATA(  I  )  ,VOLT(  I  ) 

EtiO  NEaT  1 
b  A0  FOR  f=l  TO  N 

700  PR! MTU  I  ;  OHMS_DATA(  I  )  , VOL  T ( I  ) 

7  I  0  ME  T  I 

720  ASSIGN#  1  TO  • 

730  GO  SUB  Pl.UT_RQUT  INF 
740  GOTO  770 

750  ENTER  MaME$  i  NEW  TITLES 
760  Oulu  4 til) 
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770  LUU 

vuo  n r £p  "WHnr  is  run  desired^" 

730  PEER  I  0  ,500 

610  rSSIGNtt  1  TO  "•  INSERT  NAME  OF  RUN  • " 

S20  FOR  1-1  TO  N 

POOP LmDU  1  ;  GHMS_GATA(  I  )  ,UOLT( I  ) 

3-40  NE/'l  1 
330  r.bSIGMtt  1  TO 
BE0  GOTO  ^30 
3b  I  ! 

8  62  * 

370  PLGr_RQUT JNE :  •  ***  PLOTS  POLnR  DATA  »♦* 

87  1  t 
6  72  i 

630  PLOTTER  IS  70S 
630  GRAPH 
900  Gi  l  EAR 
910  REN  i 

320  Si  hi  E  B0  ,60,-60  ,60 

930  XhA IS  0,10,-40,40 

94  0  Y  ri  *  I S  0,10,-40,40 

950  CSIZE  3 

960  LINE  1 vPE  1 

970  FOR  L*0  TO  40  STEP  10 

930  MiRE  0  ,L 

990  t  ORG  b 

1OO0  Lrf'EL  L-40 

1010  NEXT  L 

1020  l  if  if.  TYPE  1 

1030  C3I2E  6 

104  0  MO‘'E  0,50 

104  1  i 

1042  i  INSERT  LABELS  1050,  1060  AND  1100  PRIOR  TO  START  OF  RUN 

1043  I 

1050  LABEL  ; “OUADRUPOLE  BEAM  PATTERN  <dB)” 

1O60  SEIZE  3 


1 070 

M'V't  0, 

45 

1  Oj0 

1  rip L  1  ; 

“  R  -3 . 20  N "  ,"  “/I 

d -3 . 30"  , 

"  "  , "  FREQ  =  300  Hr" 

1090 

MiF'E  0, 

-GO 

1  1  OO 

l.rBFI  ; 

“ PHASE=0 ,180,0,18 0 

Deg "  , 11 

"  ,"ANPL-  0.093  Uac ''  , " 

“  , " AN6LE  =  30  Deg 

1112  NR’L  45,0 
I  1 20  LAP  0  5 
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I  13d  LABEL  ;  "MIKE" 

1140  LINE  TYPE  4 
I  I  Ed  FOR  L  = I  0  TO  40  STEP  10 
MOO  I'lOUE  I  ,0 

1170  FOR  1=0  TO  3E0  STEP  5 
I  160  DRAW  Lm'.uS  (  I  >  ,  L  »  5 1 N  (  I  ) 

I  I  90  NEXT  I 
1200  NEXT  L 
1210  LIME  TYPE  I 
1220  PEN  UP 
I230  FEN  2 

1240  FOR  1=1  TO  N  i  BEGINS  BEAM  PATTERN  COMPUTATION 
1250  IF  UOLTU)  MUuLT  THEN  MUOL  T  =  U0LT (  I  ) 

I 2S0  NEXT  I 

1270  FOR  1=1  TO  N 

1290  K=OHMS_OATA1 I  )/ 10000*300 

1290  h I =K 

1300  If-  OHMSJIATAI  I  )  1 0  THEM  K=KH-360/295 

13I0  IF  OHNS_DATA( I  )'  10500  THEN  GOTO  1400 

1220  NEC  I  0A  =  20*L6T  (  UOl.  T<  1  >/MUOl.T  )U0 

1330  LINE  TYPE  3 

1343  IF  UECTOR  I  THEN  UECTQR=0 

1350  PRINT  I  ,l-  , UECTOR 

I  3E0  Xi  I  )=*UECTOR *COS  ( K  ) 

1370  V  (  I  )  =  UEC  T  OR  *  S  IN  <K> 

1 380  PLOT  '< (  I  >  ,Y(  r  ) 

1390  PEN  UP 
1 400  NEXT  I 
1410  RETURN 
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TABLE  III 


SAMPLE  PROGRAM  DATA  COLLECTION  OUTPUT 


3 

4 

c 


OHIiSJnTA  (OHfiS) 
4 


3G 

56 

79 


MOLT  (M5C; 
.  54224 
. 70029 
.75352 
.70253 
.55367 


6 

7 

8 

9 

10 
1  1 
12 

13 

14 


106 

.61526 

123 

.57768 

1  60 

.54098 

1  94 

.52056 

215 

.51957 

251 

.53613 

281 

. 567 

306 

.60786 

330 

.65774 

369 

.72735 

16 

17 

18 
5  9 
20 

2  1 
n  -t 

n  ~7 

24 

25 

26 
27 
26 

29 

30 


34 


401 

. 50624 

429 

.66951 

456 

.93089 

430 

1  .00539 

522 

1.07149 

543 

1  .  1 2759 

572 

1  .  169^7 

600 

5.244 

623 

1  .20525 

6  4  6 

1 . 3225 

cn 

L>J 

1.36172 

704 

1 .40042 

726 

1 . 42022 

755 

1  . 4  SB  4 

706 

1  .47033 

814 

1  .49626 

646 

1 . 50758 

06  9 

1.51197 

095 

1.511 00 

924 

I  .50402 
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l.j  l-j  oj  r  j 


TABLE  IV 


SAMPLE  CAinJLATED  BEAM  PATTERN  OUTPUT 


I 

Degrees 

Decibels 

T 

1  .  Jb433S9830S 

3  1  . 3  I  7  j 3378 2 7 

J 

.  62  3 

34 . 25425^3403 

Sj 

1  .296 

33 . 954  4  396329 

4 

2.015 

33.345 1 423354 

5 

2.844 

32.5507973377 

6 

3.615 

32 . 2325595096 

7 

4.503 

31 .6428331572 

6 

5. 75 

31 .0727607057 

9 

5.354 

30.739552392 

10 

7.74 

30.7220173563 

1  1 

9.035 

30 . 9945337033 

1  2 

10.115 

3 i . 4  56004  604  3 

1  3 

11.063 

32.0852078323 

1  4 

11.33 

32.770221597 1 

15 

13.284 

33.6440053909 

1  6 

14.436 

3  4 .5354233377 

1  7 

15.444 

35. 1946381359 

1  6 

16.438 

35.7371035167 

1  3 

17.556 

36.4555277425 

20 

1 8.792 

37. 00559596 1 6 

2  1 

19.34 5 

Jj  i  .  4  b  2  1  b  tj  w  J  1  J 

20.592 

37.9163367977 

2  3 
j 
0 


34 

35 


2  1.5 
22.503 
23.323 
24  .223 


35.305544124b 
33 .5091 395402 
38 . 837050  1  1  5 
39.09059264 


25.344 

39.3343026072 

25.136 

39.5050387255 

27.13 

33.67  4  7  4  99..Z7 

28.296 

39. 504564321 6 

29.304 

39.9092773356 

30.456 

39.3747436543 

31  .284 

40 

T  “5  TO 

39. 5943d 55 571 

.  j  j . 26  4 

39.9558276 1 
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RIPOLE  BEAM  PATTERN  (dB) 


APPENDIX  C 


Experimental  and  Theoretical  Plots 
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THETA  VS.  DIRECTIVITY 


pin=9o 


•j  vjhny 


kD1=kD2=2kD3-  5.276 

Dl=0.3  D£=03 
D3-0.15  D4=0. 
P2=1B0  P3=  0  F4=  0 
j  A=l.  B=l.  0=0. 


Figure  C.2  Theoretical  Three  Element  Isosceles  Plot 
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QUADRUPCLE  BEAM  PATTERN  (dB) 
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THETA  VS.  DIRECTIVITY 


PHI=  0  DEGREE 

kDl=kD2=kD3-kD4=  5.276 

Dl=0.5  DC-05 
D3=05  D4=05 
P2=  0  P3-100  P4=180 
A=l.  B*l.  C=l. 


Figure  C.4  Theoretical  Planar  Quadrupole  Plot 
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TRIPOLE  BEAM  PATTERN  (dB) 
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TRIPOLE  BEAM  PATTERN  (dB) 
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TRIPOLE  BEAM  PATTERN  (dB) 
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QUADRIPOLE  BEAM  PATTERN  (dB) 
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GUADRUPOLE  BEAM  PATTERN  (dB) 
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